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Abstract 

Little Higgs models provide a natural explanation of the little hierarchy between 
the electroweak scale and a few TeV scale, where new physics is expected. Under 
the same inspiring naturalness arguments, this work completes a previous study on 
lepton flavor-changing processes in the Littlest Higgs model with T-parity exploring 
the channel that will eventually turn out to be the most sensitive, ^ — e conversion 
in nuclei. All one-loop contributions are carefully taken into account, results for the 
most relevant nuclei are provided and a discussion of the influence of the quark mixing 
is included. The results for the Ti nucleus are in good agreement with earlier work by 
Blanke et al., where a degenerate mirror quark sector was assumed. The conclusion 
is that, although this particular model reduces the tension with electroweak precision 
tests, if the restrictions on the parameter space derived from lepton flavor violation 
are taken seriously, the degree of fine tuning necessary to meet these constraints also 
disfavors this model. 
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1 Introduction 



Lepton flavor is conserved within experimental limits in all known processes, except for 
neutrino oscillations. This is predicted by the Standard Model (SM) when it is minimally 
extended to include Dirac and/or Majorana masses [1]. In such a case any neutral lepton 
flavor transition is proportional to the neutrino masses rrii, ~ 0.1 eV, and then negligible if 
the relevant energy scale in the process M is near or above the GeV. Hence, the observation 
of lepton flavor violation (LFV) in processes like /i — ?■ e7, /i — ?■ eee or /i N — e N, or the 
analogous r decays would be a clear signal of new physics beyond the SM. 

On the other hand, any new physics near the TeV scale involving leptons may induce 
large LFV transitions if it is not aligned with the SM charged lepton Yukawa couplings. 
In this case the relevant scales in the process are the electroweak scale v ~ 246 GeV and 
the new scale / of few TeV, with the corresponding ratio v/f ~ 0.1 almost nine orders of 
magnitude larger than rrii, /mE Due to very stringent bounds on LFV processes (as small 
as 10~^^ for n — )■ e7) new physics at the TeV must incorporate a very efficient lepton flavor 
suppression mechanism to agree with current and foreseen limits on LFV transitions. 

Present limits on flavor violating processes set stringent constraints on the possible 
extensions of the SM, in particular on those solving the (little) hierarchy problem like 
supersymmetry. Little Higgs models and models with large extra dimensions. (See for a 
recent review [2], and references therein.) These bounds, which in the quark sector and for 
r leptons have been improved by BaBar and Belle [3], are also expected to be improved 
for the first two lepton families in the near future. Thus MEG will improve the precision 
on /i — )■ e7 and fi — >■ eee by two orders of magnitude [HEj, whereas PRISM/PRIME could 
reduce the bounds on /iN — )■ eN by several orders of magnitude [6]. Super-B factories 
under study should also improve the precision on r decays by an order of magnitude [7] . 

Little Higgs models take care of the large top corrections to the Higgs mass by making 
the Higgs a pseudo-Goldstone boson of a larger global symmetry broken at a scale of a few 
TeV P4T3]. This appealing solution to the little hierarchy problem encounters difficulties 
when a specific model is defined. In particular, electroweak precision data [T3HT9] and the 
limits for fiavor violating processes tend to banish the new physics scale to higher energies. 
To deal with the electroweak data constraints, T-parity [201422] was introduced with the 
Littlest Higgs model with T-parity (LHT) being the most popular solution for a lighter 
new scale. 

In this paper we re-evaluate /iN — eN in the Littlest Higgs model with T-parity [23]. 
This complements our previous calculations of // — >■ e7 and ^ — > eee in [21] and those 
in [23l|25]. The results for the Ti nucleus are in good agreement with previous work [23|[26] 
where a degenerate mirror quark sector was assumed. Here we have also considered the 
effect of a more general quark sector where we include non-degenerate mirror quark masses 

* Note that amplitudes must be at least suppressed by ouc of these ratios, and then cross sections are 
at least quadratically suppressed. 
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and arbitrary mixings. We have also found that /iN — )■ eN is one-loop finite [2 ^ 126 1 [27]. 
This is apparent for box diagrams but not for Z penguins. In general, this process provides 
a stronger bound on the new physics scale /, but limits from all three processes are required 
because there are accidental suppressions of any of them in certain regions of parameter 
space. Obviously, each of them cancel when the heavy lepton Yukawa couplings align with 
those of the charged leptons, that is the amplitudes vanish when the corresponding mixing 
angle Q goes to zero. The amplitudes also scale like when the new physics scale goes 
to infinity with all other parameters finite. Present bounds on these processes are satisfied 
for Q < 0.01 or / > 10 TeV. Particular cancellations result from amplitude contributions 
of different sign, for example between penguins and boxes, at given parameter values. As 
we shall discuss, in the case of N — )■ e N these cancellations occur even for degenerate 
quark masses (or just one quark family). 

On the other hand, box as well as 7 penguin contributions are at most logarithmically 
divergent when the internal fermion masses go to infinity, whereas the Z penguins add a 
quadratically divergent contribution when these masses are parametrically increased. This 
is similar to what happens in the SM with the top mass lUt [281129] (and to /x — >■ eee in 
these models [21] )• In our case the heavy fermion masses are proportional to k/, where 
K is the corresponding Yukawa coupling. For a fixed / value, the Z penguins scale like 
K^, the same as, for instance, the top quark contribution to Z\)b that scales like A^, with 
rut ~ \v [30]. In both cases the gauge symmetry, [SU(2) x U(l)]i x [SU(2) x U(l)]2 in 
the LHT and SU(2)j;^ x U(l)y in the SM, is spontaneously broken and the heavy fermion 
mass splits the fermion multiplet. Of course, these Yukawa couplings can be large and the 
corresponding contributions dominant, but they should eventually remain perturbative. 

In a companion paper we present the evaluation of the three processes \i — >■ e7, \i — >■ eee 
and /iN — )■ eN in the Simplest Little Higgs model [3l], finding similar results [32]. These 
calculations can also be applied to the four-dimensional formulation of extended models 
with large extra dimensions [33l[3l]. In this case, however, care must be taken to include 
the appropriate number of Kaluza-Klein modes to match the five- dimensional result [35] . 

In Section 2 the LHT model is briefiy described with the purpose of introducing notation 
and conventions, providing the Feynman rules for quarks, needed to perform the calculation 
of /iN — J- eN which were not included in [21]. (They agree with those in [SB] up to sign 
conventions.) In Section 3 we detail the different contributions, finding complete agreement 
with [23] in the total sum. We show in Section 4 that, in general, the conversion rates for 
both Ti and Au give a more stringent limit than /i — )■ e7 and /i — )■ eee, with Au being 
the most stringent. Results for Pb are not shown but it was found that it provides a 
lower conversion rate than Ti. We also study the effect of the heavy quark sector on the 
conversion rate. A degenerate heavy quark sector was previously assumed in the literature 
arguing that the quark mixing is a higher order effect. We analyze the degenerate case, 
and then show the impact of a non-degenerate heavy quark sector including quark mixings. 
Even assuming degeneracy, the quark mass parameter can entirely cancel the conversion 
rate if chosen appropriately (in a region where box diagrams cancel penguin diagrams). 
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Although the effect of the mixings is found to be less important than that of the quark 
masses, they can have sizable effects on the conversion rate since they can push the values 
of the form factors into regions with cancellations. Finally, Section 5 is devoted to our 
conclusions. 



2 The LHT model 



The Littlest Higgs [13] is a non-linear a model based on the spontaneous breaking of a 
global SU(5) symmetry into SO (5) by the vacuum expectation value of a five-dimensional 
tensor field. The 14 Goldstone fields are parameterized as 
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Only an [SU(2) x U(l)]i x [SU(2) x U(l)]2 subgroup of SU(5) is gauged, with generators 
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and 0"" the three Pauli matrices. This gauge group is broken by Eg down to the SM 
group, whose generators are the combinations {Q1 + Q21 ^1 + ^2}- A discrete T-symmetry 
is introduced that exchanges the gauge fields Gi^2 of [SU(2) x U(l)]i^2 under the 

assumption that g = gi = g2, g' = g'l = g2- The T-even combinations remain massless 
while those T-odd acquire a mass proportional to /. In order that the SM Higgs doublet 
contained in 11 is T-even and the remaining Goldstone fields T-odd, the T transformation 
on the scalar fields is defined as follows, 



n 



-nun, n = diag( 



-1,-1,1, 



■1) 



and then 
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The gauge- and T-invariant Lagrangians for the vector and Goldstone bosons can then 
be written in terms of covariant derivatives. Our sign conventions are those in [23]. After 
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the electroweak symmetry breaking the photon remains massless, the weak gauge bosons 
pick up a mass of order the T-odd gauge bosons get rotated to their physical states 
Ah, Zh and W^, and the scalar fields are also shifted and rotated accordingly. The latter 
include the would-be Goldstone fields r], u^, eaten by the four heavy gauge bosons 
above. 

The introduction of fermions in the model is less straightforward. Since the imple- 
mentation of the lepton sector was described in detail in [21], which was sufficient for 
the processes yU — )■ e7 and fi — eee, we present below the relevant Yukawa and gauge 
interaction Lagrangians for the quark sector, also needed to study /i N — )■ e N. 



2.1 Quark Lagrangians 



As for the lepton sector, one introduces two left-handed fermion doublets in incomplete 
SU(5) multiplets, one transforming only under SU(2)i and the other under SU(2)2, for 
each SM left-handed quark doublet [22] : 
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. An SU(5) transformation V and a T transformation acts as follows: 
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The T-even combination \l/i+r2So\l/2 remains massless while the T-odd ^ i—VLT^q^ 2 obtains 
a mass when coupled (see below) to a multiplet ^ r defined as 
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so that R transforms as a 5 representation. The operator f/ is a non- linear transforma- 
tion depending on V and the Goldstone fields, that takes values in the Lie algebra of the 
unbroken SO (5), defined by 



0; 



such that the following Yukawa Lagrangian is SU(5) and T invariant, 

= -f^f {^2^^R + ^iSoe^^i?) + h.C. . 



(2.9) 
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Defining the T-eigenstates 

qiL - q2L qiL + q2L , /o 11 N 

qL = — , qHL = — , q = u,d, (2.11) 

we see tliat tlie left-handed T-odd fields uhl and duL couple to uhr and duR, respectively, 
to get masses of order /. The T-even combinations ul and di are the SM left-handed 
quarks. The fields xr and ifjR can be assumed to acquire large Dirac masses pairing with 
additional fermions [22|[37] so that they decouple from the theory. 

The left-handed quark fields interact with the gauge bosons through the following 
covariant derivatives: 

Cfl = i'l'i^i^i + i^'2^2*2, (2.12) 

where 

D^^ = d, + V2ig{WlQf + W^^Qf) + V2ig\yt^B^^ + y^^B2^), (2.13) 
D2, = d,-V2ig{W^^Ql + W^^Q^,) + V2ig'{yf^B,^ + y^^^B,,), (2.14) 

with ly/^, Bj^ the gauge fields of [SU(2) x U(l)]j and 

y?' = y!' = j^. y!' = y?' = ^- (2.15) 

It is necessary to enlarge SU(5) with two extra U(l) factors to introduce these hypercharges 
[27] . that add up to 1/6 (the hypercharge of the SM left-handed quark doublet). 

The gauge interactions of the heavy and light right-handed fields are not necessary for 
our calculation since they turn out to be of higher order in v/f, but we include them for 
reference 



Cfr = i-^R-f'' 



"^R, (2.16) 



2 ' ^ ' 2 

J^'fr = mRl''{d^ + ig'yuB^)uR + \dR-i''{d^, + ig'ydB^,)dR, (2.17) 

where yu = 2/3, yd = ~l/3 are hypercharges that also require the two additional U(l) 
factors [27] . 

In the case of three generations, the Yukawa Lagrangian ( ]2.10p that provides the masses 
for the heavy quark states is given by: 

Cy„ = -Kmnf in^n + ^T^O^^n) + h-C, (2.18) 

with family indices m^n = 1, 2, 3. This is the primary source of heavy flavor mixing in the 
quark sector, which is completely analogous to the lepton sector. After diagonalization. 
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diag(Km) = V^^kV^ where {^11) acts on the left (right) handed fields, one finds that 
the masses for the heavy quarks are approximately given by: 



Then the T-odd gauge boson interactions from f l2.12p work out as: 



(2.19) 



+ h.c., 



(2.20) 



where Gh is a T-odd gauge boson {Ah, Zh, Wh) and the matrices V^'*^ are rotations 
necessary to make the light sector mass diagonal. The matrices appearing in the coupling 
can therefore be defined as [391: 



Vhu = V^^, VHd = v^'vt 



(2.21) 



The two light rotations are then related by V^'^V^ = Vckm and can be obtained from the 
Yukawa Lagrangians for the light sector which we now describe. 

For the SM quark Yukawa couplings, the first two generations are treated separately 
from the third one. This is because of the special structure necessary for the collective 
symmetry breaking mechanism. The third generation requires enlarging and \l/2 into 
multiplets of the SU(3)i and SU(3)2 subgroups of SU(5) [39] adding additional T-even and 
T-odd partners (even and odd combinations of t[ and ^2); usually denoted by T+ and T_ 
respectively. 



Ql = (i<,-i<,t;^,o,o)'^, Qi = (o,o,4^,i4^,-iMy^. 
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They have the same transformation properties as \l/i_2 in (12. 7p . The quarks masses of this 
generation originate from the Lagrangian 
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where k = 1,2,3 and x,y = 4, 5. 

The first two generations are implemented differently and do not require the enlarged 
SU(3) multiplets: 
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■J=f^ijk^xy <5liSja;Sfcy — ((52*So) j Sj^Sfcy + h.c. (m— 1,2), (2.24) 



where k = 1,2,3; x,y = 4,5 and 



Qr = (idr^,-i<^, 0,0,0)' 



Q^ = (0,0,0,ici™^,-in^J 



(m = 1,2). 



(2.25) 
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Unlike the up quark case, all three generations of down-type light quarks acquire mass 
through the same type of Lagrangian. For simplicity, we have introduced the quark mixing 
(CKM matrix) in this sector only, choosing the up-type quarks to be directly in the mass 
basis. This is because of the fact that the three generations of up quarks are implemented 
differently, which explicitly breaks flavor symmetries. Although it is possible to restore 
flavor symmetries (see [39] for further details), this is unimportant in our case and it is 
sufficient to assume that all the mixing comes from the rotation in the down sector. Then 

d^ + h.c. (m,n = 1,2,3), 

(2.26) 

where i,j = 1, 2; x,y, z = 3, 4, 5, X = (£33)^^/'^ and 

VI/7 = (n- , rf™ , 0, 0, 0)^ , = (0, 0, 0, d^^f . (2.27) 

In this context, the previously deflned is set to the identity matrix since there is no 
family mixing in the up sector and would correspond to Vckm and would be one of the 
two matrices necessary to diagonalize A™". The other matrix, would be the rotation 
for the right handed sector. 

The relevant degrees of freedom in Vhu (or Vud) are three angles, 9^2, ^23 and 6*5*3, 
three complex phases, as emphasized in [38] where a standard parameterization is proposed. 

2.2 Feynman rules for the quark sector 

From these Lagrangians one can obtain the Feynman rules for the quark vertices entering 
in our calculation. They are in agreement with [36] and are summarized in tables [1] and [21 
in our sign conventions. Standard Model Feynman rules are used for vertices that involve 
only ordinary particles and they will not be listed. Although the model predicts some 
corrections to some of these standard vertices, their contributions are subleading in our 
process and we can safely ignore them. 

3 Contributions to the /i N — > e N process 
3.1 Topologies and form factors 

The N — )■ e N process receives contributions from both penguin and box diagrams like 
those shown in flgure [Dd. Here the quark q is either an up or down quark and is of the 
same type before and after the interaction. The gauge bosons exchanged in the penguin 
diagrams are a Z boson or a photon. The diagrams look very similar to those involved in 
fi — )■ eee except for the fact that now up and down-type quarks instead of electrons are 
involved on the lower legs. 
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Table 1: VFF vertices 'ie'~^^{gLPL + QrPr) for quarks in the LHT model; xh = 
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Table 2: SFF vertices ic^clPl + crPr) for quarks in the LHT model; xh 
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-I q 




(a) 

Figure 1: Vertex subdiagrams (a) and topologies (b) involved in /xN — )■ eN. 

The contributions of the lepton flavor- changing vertex subdiagrams fit into the following 
Lorentz structure [lDH12]l!l 

iT^{puP2) = ie [YiFlPL + FlPn) - [iFl + FW)a>^''Q, + (iFj + Fl^'^)Q^] , (3.1) 

where V denotes the external vector boson, either Z or 7 in our case, and Q = Pi — P2 
as in figure [1^. The penguin diagrams for /i — )■ eqq can be read from [23] introducing the 
corresponding electric charges {Qg) and weak couplings {Zl j^) for quarks: 

A^7Peng = [Q^riA^PL + A^Pr) + ima^'' Q .{A^ Pl + A^Pr)] u{p) 

XMg(p2)Qg7M^g(P3), (3.2) 
-Mzpcng = ^M(pi)[7^(FiPi + FRP^)]M(p)M,(p2)[7M(^l^L + ^^^R)]^'.(P3), (3.3) 

z 

where m is the muon mass, the electron mass has been neglected, and 



A{ = F2/Q', = F^/Q^ = + iF^)/m, A^ = -{Fl - iF^)/ 



Fl = -FI, Fr = -FI (3.4) 

The form factors Fs and Fp do not contribute in any case since they are found to be 
proportional to the ratio of the external lepton and quark masses to the mass of the W 
boson and are therefore negligible. The same occurs in the case of Fe and Fm in the Z 
boson penguin. 

Similarly, in the limit of vanishing external momenta the contribution to the amplitude 



T A minus sign has been extracted from the dipolc form factors Fj^i and Fe in order to meet the 
definition of the momentum transfer Q in |40j , that had opposite sign in (4HI42| . This implies a sign flip 
in ()3.4|) accordingly. 
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Nucleus 


Z N F{q) r,,pt [GeV] 


iTi 


22 26 17.6 0.54 1.7 x lO'^® 
79 118 33.5 0.16 8.6 x lO'^^ 



Table 3: Relevant input parameters for the nuclei under study. From 

from any box diagram can be written as |40j : 

Ml^ = e'B^^ [uipih'-PMp)] Mp^h.PLvM] 

+e^B^g [u{pi)YPlu{p)] [uMl^PRVgip^)] 

+e^B^^ [u{pi)j^'Pnu{p)] [uMi^PlvM] 

+e^B^^ [u{p^)PLu{p)] [ug{p2)PLVq{p3)] 

+e^Bfg [u{p^)Pnu{p)] [uMPrvM] 

+e^B^g [uip^)a''•'PMp)] [uMa^.PLVgips)] 

+e^B^^ [uip,)a^''PRuip)] [uMct^.PrvM] . (3.5) 

However, we may restrict ourselves to the first term, proportional to B^^, since it is the 
only contributing due to the fact that the LHT couplings are primarily left-handed. 

Finally, in the process /iN — >■ eN only the quark vector current contributes [l3]. The 
conversion width is then given by [ID]: 

r(;, ^ e) = a'^Fiq)'ml \2Z{A{ - A^) - {2Z + N)BI - {Z + 2iV)5f/ , (3.6) 



where 



B',, = + ^^{Zl + ZD ^ B{^ + Fl^ + Fl^. (3.7) 
The conversion rate is obtained by dividing by the capture rate: 

n = ^^. (3.8) 

capt 

The nuclei we will consider are g^Ti and ^ygAu, whose relevant parameters are listed in 
table El We now calculate the form factors appearing in the amplitudes in terms of standard 
loop functions. 



3.2 Penguin diagrams 

For the penguin diagrams, the contributions to the form factors are identical to the /i — eee 
case since they only depend on the lepton triangle diagrams. The only difference here comes 
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from the couphngs to quarks and these are taken into account in the expressions (13 ■2p and 
(13. 3 p for the matrix elements. We therefore just recall our previous results^ 



aw 1 v"^ 
aw 1 v"^ 



aw 



Stt swCw 8/^ 



Y.V-V;,^,[G^Siy^) + Gf{y^) + \Gf{ayS\, (3.9) 



E ^^H'*^^Hi [Fw{y^) + Fz{yi) + \Fz{ay,)\ , (3.10) 



Y,Vr,Vil',y,Hw{y^). (3.1i; 



where 



Hi 



a 



w 



Ml 



2 ' 
^W 



(3.12) 



and the loop functions are given in equations (3.23), (3.29), (3.35), (3.39) and (3.42) of [2^ 



3.3 Box diagrams 

The boxes do differ considerably and need to be recalculated. The list of diagrams is shown 
in figure [21 The contributions are split into two form factors depending on the quark that 
enters on the lower legs: B^^ and B^^. Notice that all diagrams with neutral gauge bosons 
include several contributions, i.e. diagrams with one or other type of boson and diagrams 
with mixed types. The first diagram in the second row, for instance, includes four of them: 

The resulting form factors follows: 



B 



aw 1 1 V \ - „ 

''J 



4 + -y,y] doiyi, y]) + 2yiy]do{yi, y 



3 ~ 3 ~ 3 ~ 



(3.13) 



aw 111' V ^ ^ 



-Aoivuyf, 



yiyj doiyi, yj) - 2yiyJdo{yi, yj) 



do{ayi, ay'f) - —doia, ayi, ayf) 



100a 



10 



(3.14) 



t The sign of the dipole form factor has been reversed as compared to [23] , where the sign of Q was 
mistaken. This has no consequences in — )■ e7 and a non-significant impact on /x — )■ eee. 
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Here we have introduced 



Xij — ^HrHe l^Hul ' Xij — ^Ht^Hl \ ^Hd\ ' {6.10} 

and 

2 2 



The box functions (io(x, ?/, 2), (io(x, y, z), do{x, y) and (io(x, y) are given in equations (C.26) 
to (C-29) of 



4 Numerical Results 

Next we can proceed to compute the conversion rates for different values of the LHT 
parameters. As in our previous work [23], we restrict ourselves to the case of two lepton 
generations. Then, the lepton sector can be fully determined by four parameters: The LH 
breaking scale /, the masses of the heavy lepton^ rnui {i = 1,2), and the angle 9 that 
defines the now 2x2 mixing matrix in the lepton sector Vne'- 

V,,= l ^«^\^{ I (4.1) 

Vm Vm j \ -sinS cosO ' ^ ' 

We shall again replace the masses mm by the parameters 5 and y defined as: 

2 

y = Vym, yr = T7F^ ^ = 1^2, (4.2) 

S = ""^^'""^^ (4.3) 

Both niHi and Mw^ are proportional to / and, therefore, y and S are independent of this 
scale. 

All form factors for this process then take the following general form: 

E n7V}r,F{y,) = '-^[F{yi) - F{y,)], (4.4) 

i=l,2 

with F a generic function. The dependence of the conversion rate for small 6 can be 
approximated by: 



7^ oc 



5 sin 2^ 

J 



2 



(4.5) 



^ We neglect the difference between the masses of the heavy neutrino and the charged lepton of the 
same generation, that is of order of (u//)^, as for quarks in (|2.19p . 
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Figure 2: Box contributions to /iN — )■ eN in the LHT model. 
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Here the dependence on / and 6 is exact. There is, of course, a dependence in y but the 
behaviour with changes in this parameter cannot be expressed as simply. 

The quark sector requires another set of masses and mixings. The two mixing matrices 
involved in the quark sector, Vhu and Vnd, are related by Vnd = VhuVckm so there is, 
in fact, only one matrix that requires fixing. We will assign masses to the three heavy 
quarks ignoring the small mass difference f l2.19p between up- and down-type quarks of the 
same generation (yj* ^ yf). The parameters for this sector are defined analogously to the 
three-family lepton sector in 



2 

m i 

r = Vyfy^^ y'i = T^^ ^ = 1,2,3, (4.6) 

2 2 2 2 

m 2 — TTi 1 m 3 — m 2 



4.1 Case without quark mixing 

We will firstly assume no flavor mixing in the quark sector, which can be achieved by 
assuming mass degeneracy of heavy quarks. This is convenient to view the general behavior 
of the model clearly. In figure [3] we show the conversion rates divided by the current 
experimental limits (7^(^uAu eAu) < 7 x IQ-^^ |15] and 7^(/iTi eTi) < 4.3 x IQ-^^ 
|46]). Only values below unity are experimentally allowed. The results of [23] for fi 
{B < 1.2 X 10"^^ [17]) and fx eee {B < 10"^^ [48]) are also plotted using the same 
inputs. The masses for the heavy quarks have been chosen to correspond to = 1 and 
6^2 = ^23 — 0- Natural values for / = 1 TeV, lepton mixing 6 = tt/A, mass splitting 6 = 1 
and y = 4 are taken as a reference, here and in the following unless stated otherwise. Notice 
that the reference value for y has been shifted from the one used in [23] . The reason for this 
can be seen in the top right hand graph of figure [3] which shows the dependence on y. The 
value y = 1 chosen in [23] falls very close to a point where penguin and box diagrams cancel 
each other out. This artificially reduces the value of the conversion rate to values below 
the values expected for — )■ eee and yU — )■ e7 even though in almost any other region the 
/i N — )■ e N process dominates. These cancellations also occur elsewhere in the parameter 
space and in the other processes. This only means that the process in question does not 
set a bound on the parameters in that area. The key observation is that the rates for the 
three processes, // — )■ e7, /i — )■ eee and yU N — )■ e N, in general do not simultaneously vanish 
for any of the values of the parameters making it impossible for the model to fit into those 
areas. Scanning over parameter space has been discussed in some detail previously [23||26]. 
finding large correlations for non- vanishing cross-sections between different processes. Here 
we prefer to look for the parameter regions where all processes get small, fulfilling all 
experimental constraints. This is so when the heavy Yukawa couplings effectively align 
with those of the charged leptons, i.e. small 6 or 6. Obviously, all new contributions also 
cancel for a large new physics scale, since the corresponding cross-sections scale like f~^. 
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Figure 3: Ratios of LHT predictions to current limits, assuming = 1 and degenerate 
heavy quark masses, as functions of the different parameters keeping the others at their 
reference values. 
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/i — > e7 


fi — )■ eee 


/i Au — e Au 


/i Ti — > e Ti 


Limit 


1.2 X 10-" 


10-13 


10-12 


10-14 


7 X 10-13 


4.3 X 10-12 


10-18 


//TeV > 


3.00 


9.93 


3.98 


12.6 


9.11 


5.13 


234 


sin 26 < 


0.111 


0.010 


0.055 


0.006 


0.012 


0.038 


< 10-4 


\6\ < 


0.106 


0.010 


0.062 


0.006 


0.012 


0.039 


< 10-4 



Table 4: Constraints on LHT parameters from present and future experimental exclusion 
limits, assuming = 1 and degenerate heavy quark masses. 



Barring the possibility of these cancellations, it is clear from figure [3] that the most 
restrictive process is /x Au — )■ e Au, whose constraints are somewhat more demanding than 
those previously obtained in [24] for — e7 and fi — )■ eee and than the ones for /iTi — )■ 
eTi. However, this last process is expected to have the greatest improvements in future 
experiments |1H6]. We compare the bounds on the parameters derived from n ^ e'y and 
H — >■ eee with the new ones from /i N — >■ e N in table |4]for current and future measurements. 
The bounds on each of the parameters come from keeping all the others at the reference 
values and finding the region where the conversion rate (or branching ratio) is within the 
experimental limit. Notice that these bounds depend strongly on the choice of input values 
so these bounds are not strict. 

A comment about the large y behavior is in order. Figure H] shows the contributions of 
photon penguins, Z penguins, boxes and their coherent sum to the total /i — e conversion 
rate in Au. Notice that only the Z penguins are responsible for the non-decoupling effect, 
similar to that of the top quark in the SM [281430] . They saturate the /i — e conversion 
rate for large y, analogously as for fi — )■ eee in figure [31 For such large values the main 
contributions results from fl3.1ip . which grows like k,iK2 for large heavy Yukawa coupling 

i^i, y = y/mrn oc kiK2- 

Figure [5] shows exclusion contours in the (sin 2^, 5) plane for the present experimental 
constraints for fi Au — )■ eAu. The case of y = 1 shows large allowed areas due to the fact 
that this is close to a cancellation point. In the other cases, the mixing angle and the mass 
splitting are correlated if we are to remain within the experimental bound, as they were 
for /i — e7 and fi — )■ eee. The bound is relaxed for higher values of the scale /. 

In figure[6]we show the dependence on the heavy quark mass parameter while keeping 
all quark masses degenerate. We observe that the location of the cancellations depends 
strongly on the value of the lepton mass parameter y. Because of these regions, even a 
degenerate heavy quark sector can be tuned to suppress the lepton flavor changing effect. 
We therefore conclude that the values of the quark masses can have large effects on the 
conversion rate. The regions themselves, however, are relatively narrow concentrating on 
very definite values. 
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Figure 4: Contributions from photon penguins, Z penguins, boxes and their coherent sum 
to the total fj, — e conversion rate in Au, assuming = 1, degenerate heavy quark masses 
and reference values for the remaining parameters. 




-3-2-10 1 2-3-2-10 1 2-3-2-10 1 2 3 



5 6 5 

Figure 5: Contours of 7^(/i Au ^ e Au) = 4.3 x lO'^^ ^^lues of / = 0.5, 1, 2, 3, 4 TeV 
(from bottom up) and y = 0.25, 1, 4 (left to right), assuming y'^ = 1 and degenerate heavy 
quark masses. 
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Figure 6: Dependence of the N — )■ e N conversion rate on the heavy quark mass parameter 
for degenerate heavy quarks and several values of the lepton mass parameter y. 



4.2 Quark mixing effects 

To study the effect on the conversion rate of the mixing matrix as it compares to the 
effect of the quark masses, we must choose the quark mass parameters in an area free 
of accidental suppressions. One such choice is the quark mass parameter = 1 and the 
lepton mass parameter y = 4. With these values, in figure [7] we make a scatter plot showing 
the conversion rate as a function of the quark mass splittings. The quark mixing angles 
and phases in Vhu (YHd = VhuVckm) are uniformly distributed in their full range. The 
solid lines correspond to minimal heavy quark mixing, that is, alignment of heavy and light 
quark flavors {Vhu = I and Vnd = Vckm)- 

From the right panel of figure [7] we can conclude that the mixing angles, although 
numerically important, do not alter the result as much as the mass parameters y", 
and 523) long as we are far from a cancellation area. However, taking y = 1 the right 
panel shows that one can obtain virtually any value for the conversion rate by choosing 
the mixing angles appropriately. 



5 Conclusions 



As previous studies have demonstrated, LHT models are heavily constrained by flavor. 
In this paper we complete our study presented in [21] with the third basic LFV process, 
/xN— )-eN. We have complemented other calculations [23l[25] of lepton flavor violating 
effects for this process and compared with predictions for — )■ e7 and /i — )■ eee. 
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Figure 7: Ratio of conversion rate to current limit as a function of the heavy quark 
mass sphttings for random values of mixing angles. The solid lines correspond to 
VHd = VhuVckm = VcKM- The average quark mass is fixed to y'' = 1 and the lepton 
mass parameter to y = 4 (left) and y = 1 (right). Red (green) points are for Au (Ti). 



We have first briefly reviewed the quark sector Lagrangians for the LHT and find 
agreement with Feynman rules calculated in [36]. Generic limits on the LHT parameters 
can be found in table HI In earlier works only conversion in Ti was considered but we find 
that, currently, conversion in Au gives the most stringent limits on the LHT parameters 
in general, producing normalized conversion rates (i.e. the conversion rate divided by the 
current experimental limit) of up to an order of magnitude larger than in the Ti case and 
requiring a Little Higgs breaking scale of almost 10 TeV or fine tuning at the percent level 
of the mixing angles and mass splittings in the lepton sector, which must be correlated in 
order to fulfill current experimental limits. However, future limits on the Ti conversion 
rate will clearly be the most restrictive. 

Nevertheless, certain values of the different parameters can suppress the conversion 
rate by causing accidental cancellations among the different contributions. This means 
that, in some small regions, these processes do not actually provide any restrictions on the 
model, and only those derived from the other processes remain. It is important to note, 
however, that the cancellation regions for the various processes do not in general overlap 
and therefore do not allow for the model to survive within them, since there is always at 
least one process that is above the experimental limits. For /i N — )■ e N, the origin of these 
cancellations is typically a sign difference between the penguin contributions and boxes, 
which allows them to cancel each other. 

We have also considered the effect of a more general quark sector where we include 
non-degenerate quark masses and arbitrary mixings, whereas previous studies assumed 
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a degenerate quark sector. We find that even in the degenerate case, there is a sizable 
influence of the heavy quark mass parameter and, in some cases, it can completely cancel 
the conversion rate by pushing the prediction into one of the aforementioned suppression 
regions. We also observe that the remaining quark sector parameters, namely mass splitting 
and mixing angles, can also shift the predictions somewhat although, as may be expected, 
the effect is smaller (less than a factor 10 in the normalized conversion rate) as long as we 
are far from a suppression region. Otherwise, moving the angles and splittings can again 
accidentally reduce the conversion rate and change the value considerably or even cancel 
it completely. 
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A Comparison with previous calculations 



We present below the relation of the form factors for /i — )■ e7, /i — )■ eee and /i N — )■ e N used 
in this and our previous work [21], which follow from [KJIH^. with those introduced by 



Buras et al. 



corrected in 
Gf 
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^eMd 
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2V2' 



'odd 
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In the comparison with their results for the flnal expressions we have found agreement. 
Nevertheless, left and right hand sides of (lA.3HA.5l) differ by a common function of the 
heavy fermion mass involved in the line with no flavor change between external legs. 



Gf 



2^TT^ 



+ 



-l/j)liiyj 

(1 - y^f 



(A.6) 
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that vanishes due to the unitarity of the mixing matrices, 

Y,X^, = ^■ (A.7) 

i 
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